Bayesian molecular phylogenetics: estimation of divergence dates and hypothesis testing

With the advent of automated sequencing, sequence data are now available to help us understand the functioning of our genome, as well as its history. To date,powerful methods such as maximum likelihood have been used to estimate its mode and tempo of evolution and its branching pattern. However, these methods appear to have some limitations. The purpose of this thesis is to examine these issues in light of Bayesian modelling, taking advantage of some recent advances in Bayesian computation. Firstly, Bayesian methods to estimate divergence dates when rates of evolution vary from lineage to lineages are extended and compared. The power of the technique is demonstrated by analysing twenty-two genes sampled across the metazoans to test the Cambrian explosion hypothesis. While the molecular clock gives divergence dates at least twice as old as those indicated by the fossil records, it is shown (i) that modelling rate change gives results consistent with the fossils, (ii) that this improves dramatically the fit to the data and (iii) that these results are not dependent on the choice of a specific model of rate change.Results from this analysis support a molecular explosion of the metazoans about 600 million years (MY) ago, i.e. only some 50 MY before the morphological Cambrian explosion. Secondly, two new Bayesian tests of phylogenetic trees are developed. The first aims at selecting the correct tree, while the second constructs confidence sets of trees. Two other tests are also developed, in the frequentist framework. Based on p-values adjusted for multiple comparisons,they are built to match their Bayesian counterparts. These four new tests are compared with previous tests. Their sensitivity to model misspecification and the problem of regions is discussed. Finally, some extensions to the models examined are made to estimate divergence dates from data of multiple genes, and to detect positive selection

Dating Phylogenies with Hybrid Local Molecular Clocks

BACKGROUND: Because rates of evolution and species divergence times cannot be estimated directly from molecular data, all current dating methods require that specific assumptions be made before inferring any divergence time. These assumptions typically bear either on rates of molecular evolution (molecular clock hypothesis, local clocks models) or on both rates and times (penalized likelihood, Bayesian methods). However, most of these assumptions can affect estimated dates, oftentimes because they underestimate large amounts of rate change. PRINCIPAL FINDINGS: A significant modification to a recently proposed ad hoc rate-smoothing algorithm is described, in which local molecular clocks are automatically placed on a phylogeny. This modification makes use of hybrid approaches that borrow from recent theoretical developments in microarray data analysis. An ad hoc integration of phylogenetic uncertainty under these local clock models is also described. The performance and accuracy of the new methods are evaluated by reanalyzing three published data sets. CONCLUSIONS: It is shown that the new maximum likelihood hybrid methods can perform better than penalized likelihood and almost as well as uncorrelated Bayesian models. However, the new methods still tend to underestimate the actual amount of rate change. This work demonstrates the difficulty of estimating divergence times using local molecular clocks

The Metabolic Consequences of Hepatic AMP-Kinase Phosphorylation in Rainbow Trout

AMP-activated protein kinase (AMPK), a phylogenetically conserved serine/threonine protein kinase, is proposed to function as a “fuel gauge” to monitor cellular energy status in response to nutritional environmental variations. However, in fish, few studies have addressed the metabolic consequences related to the activation of this kinase. This study demonstrates that the rainbow trout (Oncorhynchus mykiss) possesses paralogs of the three known AMPK subunits that co-diversified, that the AMPK protein is present in the liver and in isolated hepatocytes, and it does change in response to physiological (fasting-re-feeding cycle) and pharmacological (AICAR and metformin administration and incubations) manipulations. Moreover, the phosphorylation of AMPK results in the phosphorylation of acetyl-CoA carboxylase, a main downstream target of AMPK in mammals. Other findings include changes in hepatic glycogen levels and several molecular actors involved in hepatic glucose and lipid metabolism, including mRNA transcript levels for glucokinase, glucose-6-phosphatase and fatty acid synthase both in vivo and in vitro. The fact that most results presented in this study are consistent with the recognized role of AMPK as a master regulator of energy homeostasis in living organisms supports the idea that these functions are conserved in this piscine model

Comparisons of host mitochondrial, nuclear and endosymbiont bacterial genes reveal cryptic fig wasp species and the effects of Wolbachia on host mtDNA evolution and diversity

Background Figs and fig-pollinating wasp species usually display a highly specific one-to-one association. However, more and more studies have revealed that the "one-to-one" rule has been broken. Co-pollinators have been reported, but we do not yet know how they evolve. They may evolve from insect speciation induced or facilitated by Wolbachia which can manipulate host reproduction and induce reproductive isolation. In addition, Wolbachia can affect host mitochondrial DNA evolution, because of the linkage between Wolbachia and associated mitochondrial haplotypes, and thus confound host phylogeny based on mtDNA. Previous research has shown that fig wasps have the highest incidence of Wolbachia infection in all insect taxa, and Wolbachia may have great influence on fig wasp biology. Therefore, we look forward to understanding the influence of Wolbachia on mitochondrial DNA evolution and speciation in fig wasps. Results We surveyed 76 pollinator wasp specimens from nine Ficus microcarpa trees each growing at a different location in Hainan and Fujian Provinces, China. We found that all wasps were morphologically identified as Eupristina verticillata, but diverged into three clades with 4.22-5.28% mtDNA divergence and 2.29-20.72% nuclear gene divergence. We also found very strong concordance between E. verticillata clades and Wolbachia infection status, and the predicted effects of Wolbachia on both mtDNA diversity and evolution by decreasing mitochondrial haplotypes. Conclusions Our study reveals that the pollinating wasp E. verticillata on F. microcarpa has diverged into three cryptic species, and Wolbachia may have a role in this divergence. The results also indicate that Wolbachia strains infecting E. verticillata have likely resulted in selective sweeps on host mitochondrial DNA

Rooting and dating maples (Acer) with an uncorrelated-rates molecular clock

Simulations suggest that molecular clock analyses can correctly identify the root of a tree even when the clock assumption is severely violated. Clock-based rooting of phylogenies may be particularly useful when outgroup rooting is problematic. Here, we explore relaxed-clock rooting in the Acer/Dipteronia clade of Sapindaceae, which comprises genera of highly uneven species richness and problematic mutual monophyly. Using an approach that does not presuppose rate autocorrelation between ancestral and descendant branches and hence does not require a rooted a priori topology, we analyzed data fromup to seven chloroplast loci for some 50 ingroup species. For comparison,weused midpoint and outgroup rooting and dating methods that rely on rooted input trees, namely penalized likelihood, a Bayesian autocorrelated-rates model, and a strict clock. The chloroplast sequences used here reject a single global substitution rate, and the assumption of autocorrelated rates was also rejected. The root was placed between Acer and Dipteronia by all three rooting methods, albeit with low statistical support. Analyses of Acer diversification with a lineage-through-time plot and different survival models, although sensitive to missing data, suggest a gradual decrease in the average diversification rate. The nine North American species of Acer diverged from their nearest relatives at widely different times: eastern American Acer diverged in the Oligocene and Late Miocene; western American species in the Late Eocene and Mid Miocene; and the Acer core clade, including A. saccharum, dates to the Miocene. Recent diversification in North America is strikingly rare compared to diversification in eastern Asia

The nonadaptive nature of the H1N1 2009 Swine Flu pandemic contrasts with the adaptive facilitation of transmission to a new host

<p>Abstract</p> <p>Background</p> <p>The emergence of the 2009 H1N1 Influenza pandemic followed a multiple reassortment event from viruses originally circulating in swines and humans, but the adaptive nature of this emergence is poorly understood.</p> <p>Results</p> <p>Here we base our analysis on 1180 complete genomes of H1N1 viruses sampled in North America between 2000 and 2010 in swine and human hosts. We show that while transmission to a human host might require an adaptive phase in the HA and NA antigens, the emergence of the 2009 pandemic was essentially nonadaptive. A more detailed analysis of the NA protein shows that the 2009 pandemic sequence is characterized by novel epitopes and by a particular substitution in loop 150, which is responsible for a nonadaptive structural change tightly associated with the emergence of the pandemic.</p> <p>Conclusions</p> <p>Because this substitution was not present in the 1918 H1N1 pandemic virus, we posit that the emergence of pandemics is due to epistatic interactions between sites distributed over different segments. Altogether, our results are consistent with population dynamics models that highlight the epistatic and nonadaptive rise of novel epitopes in viral populations, followed by their demise when the resulting virus is too virulent.</p

Episodic Evolution and Adaptation of Chloroplast Genomes in Ancestral Grasses

It has been suggested that the chloroplast genomes of the grass family, Poaceae, have undergone an elevated evolutionary rate compared to most other angiosperms, yet the details of this phenomenon have remained obscure. To know how the rate change occurred during evolution, estimation of the time-scale with reliable calibrations is needed. The recent finding of 65 Ma grass phytoliths in Cretaceous dinosaur coprolites places the diversification of the grasses to the Cretaceous period, and provides a reliable calibration in studying the tempo and mode of grass chloroplast evolution.By using chloroplast genome data from angiosperms and by taking account of new paleontological evidence, we now show that episodic rate acceleration both in terms of non-synonymous and synonymous substitutions occurred in the common ancestral branch of the core Poaceae (a group formed by rice, wheat, maize, and their allies) accompanied by adaptive evolution in several chloroplast proteins, while the rate reverted to the slow rate typical of most monocot species in the terminal branches.Our finding of episodic rate acceleration in the ancestral grasses accompanied by adaptive molecular evolution has a profound bearing on the evolution of grasses, which form a highly successful group of plants. The widely used model for estimating divergence times was based on the assumption of correlated rates between ancestral and descendant lineages. However, the assumption is proved to be inadequate in approximating the episodic rate acceleration in the ancestral grasses, and the assumption of independent rates is more appropriate. This finding has implications for studies of molecular evolutionary rates and time-scale of evolution in other groups of organisms

Plasmodium knowlesi: Reservoir Hosts and Tracking the Emergence in Humans and Macaques

Plasmodium knowlesi, a malaria parasite originally thought to be restricted to macaques in Southeast Asia, has recently been recognized as a significant cause of human malaria. Unlike the benign and morphologically similar P. malariae, these parasites can lead to fatal infections. Malaria parasites, including P. knowlesi, have not yet been detected in macaques of the Kapit Division of Malaysian Borneo, where the majority of human knowlesi malaria cases have been reported. In order to extend our understanding of the epidemiology and evolutionary history of P. knowlesi, we examined 108 wild macaques for malaria parasites and sequenced the circumsporozoite protein (csp) gene and mitochondrial (mt) DNA of P. knowlesi isolates derived from macaques and humans. We detected five species of Plasmodium (P. knowlesi, P. inui, P. cynomolgi, P. fieldi and P. coatneyi) in the long-tailed and pig-tailed macaques, and an extremely high prevalence of P. inui and P. knowlesi. Macaques had a higher number of P. knowlesi genotypes per infection than humans, and some diverse alleles of the P. knowlesi csp gene and certain mtDNA haplotypes were shared between both hosts. Analyses of DNA sequence data indicate that there are no mtDNA lineages associated exclusively with either host. Furthermore, our analyses of the mtDNA data reveal that P. knowlesi is derived from an ancestral parasite population that existed prior to human settlement in Southeast Asia, and underwent significant population expansion approximately 30,000–40,000 years ago. Our results indicate that human infections with P. knowlesi are not newly emergent in Southeast Asia and that knowlesi malaria is primarily a zoonosis with wild macaques as the reservoir hosts. However, ongoing ecological changes resulting from deforestation, with an associated increase in the human population, could enable this pathogenic species of Plasmodium to switch to humans as the preferred host

Novel transposable elements from Anopheles gambiae

<p>Abstract</p> <p>Background</p> <p>Transposable elements (TEs) are DNA sequences, present in the genome of most eukaryotic organisms that hold the key characteristic of being able to mobilize and increase their copy number within chromosomes. These elements are important for eukaryotic genome structure and evolution and lately have been considered as potential drivers for introducing transgenes into pathogen-transmitting insects as a means to control vector-borne diseases. The aim of this work was to catalog the diversity and abundance of TEs within the <it>Anopheles gambiae </it>genome using the PILER tool and to consolidate a database in the form of a hyperlinked spreadsheet containing detailed and readily available information about the TEs present in the genome of <it>An. gambiae</it>.</p> <p>Results</p> <p>Here we present the spreadsheet named AnoTExcel that constitutes a database with detailed information on most of the repetitive elements present in the genome of the mosquito. Despite previous work on this topic, our approach permitted the identification and characterization both of previously described and novel TEs that are further described in detailed.</p> <p>Conclusions</p> <p>Identification and characterization of TEs in a given genome is important as a way to understand the diversity and evolution of the whole set of TEs present in a given species. This work contributes to a better understanding of the landscape of TEs present in the mosquito genome. It also presents a novel platform for the identification, analysis, and characterization of TEs on sequenced genomes.</p